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The low velocity part of a kinetic equilibrium dark matter 
distribution has higher phase space density and is more eas- 
ily incorporated in formation of a low mass galaxy than the 
high velocity part. For relativistically decoupling fermions 
(bosons), this explains one (two) orders of magnitude of the 
observed trend, that phase space densities in dark matter halo 
cores are highest in the smallest systems, and loosens con- 
straints on particle masses significantly. For non-relativistic 
decoupling and/or finite chemical potentials even larger ef- 
fects may occur. It is therefore premature to dismiss dissipa- 
tionless particle distributions as dark matter on the basis of 
phase space arguments. 
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It has recently become clear, that the otherwise suc- 
cessful cold dark matter model (CDM) for cosmic struc- 
ture formation has several severe problems, such as pre- 
dictions of cusps in the central density profiles of galax- 
ies, too many low-mass subclumbs within dark matter 
halos, and a lack of angular momentum in galaxy disks 
compared to observations. In one way or the other these 
problems are all related to the fact that CDM has no ini- 
tial velocity spread, and therefore infinite initial density 
in phase space. Selfinteractions among the dark matter 
particles has been suggested as a possible solution [Q. 
Another suggested solution has been the reintroduction 
of warm dark matter (WDM) consisting of particles with 
a moderate primordial velocity spread Contrary to 
CDM thermal WDM particles lead to finite core den- 
sity, but apparently require a high mass and therefore 
extremely early decoupling from the primordial plasma 
to account for the observed core phase space density 
of lO~ 4 M0/pc 3 (km/s) -3 in dwarf spheroidal galaxies. 
There may also be difficulties explaining the decrease in 
phase space density by a factor 10-100 for dwarf spirals 
and low surface brightness galaxies, and a further factor 
of 10-100 for normal spirals |j||§]. A single particle mass 
for WDM would appear to lead to a definite prediction 
for the central halo phase space density because of con- 
servation of fine-grained phase space density (Liouvilles 
theorem). 

However, whereas such arguments are correct in terms 
of the average phase space density, they do not take 
into account, that while originally almost uniform in real 
space, the fine-grained density is an (exponentially) vary- 
ing function of position in momentum space. Low mo- 
mentum particles are in denser parts of phase space than 
high momentum particles, and depending on the actual 
distribution function (fermion or boson, zero or nonzero 



chemical potential, relativistic or nonrelativistic decou- 
pling), the densest part of the distribution may have 
phase space density significantly above average. Further- 
more, halo formation should typically include particles 
from the low momentum end first. A system with low 
gravitational potential like a dwarf spheroidal will effec- 
tively probe only the densest part of the phase space dis- 
tribution, whereas a large spiral galaxy probably contains 
something close to the average phase density (possibly 
diluted by mergers etc). 

Other effects (baryons, mergers, phase space dilu- 
tion during gravitational collapse, particle selfinterac- 
tions etc) can further enhance the diversity in observed 
phase space densities in galaxy cores, and as demon- 
strated in a merging hierarchy where larger systems 
are gradually formed by merging of smaller units can 
explain many of the features observed. But as demon- 
strated in the following, even in the absence of mergers, 
a significant spread in primordial phase space densities 
in dark matter cores is predicted, with core densities de- 
creasing with increasing escape velocity of the system in 
question. The magnitude of the effect ranges from one 
order of magnitude for relativistically decoupling, non- 
degenerate fermions, to several orders of magnitude for 
non-relativistically decoupling bosons. Ad hoc additions 
to these simplest distribution functions would allow a 
further range. Thus, there may still be room for a single 
dissipationless elementary particle explanation of dark 
matter in dwarf spheroidal as well as dwarf spiral galax- 
ies ||. Furthermore, the natural selection of the highest 
phase density particles in the smallest systems leads to a 
significant reduction in the minimum mass for the parti- 
cle responsible, loosening the rather strong requirements 
on the epoch of dark matter decoupling in the early Uni- 
verse. 

An isotropic gas of particles in kinetic equilibrium has 
a spatial number density 



n = {g/h 3 ) / 4Trf(p)p 2 dp, 



(1) 



where h is Planck's constant, g is the number of helicity 
states, and 



f(p) = {exp [(E - n)/hT] ± 1}- 



(2) 



is a Fermi-Dirac (+) or Bose-Einstein (— ) distribution 
with chemical potential /i. Energy E is related to mo- 
mentum p and particle mass m via E 2 = (pc) 2 + m 2 c 4 , 
where c is the speed of light. 

A particle species that decouples from the remaining 
plasma in the early Universe at temperature To red- 
shifts its momenta in proportion to the expansion of the 
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Universe, p = pe>Rd/R, where R(t) is the cosmic scale 
factor, and its number density (particle number is con- 
served) evolves like n oc R~ 3 . From this follows that the 
distribution function / at a time after decoupling is re- 
lated to the distribution at decoupling fp from Eq. (^) 
like 



f(p) = f D ( P R/R 



D) 



(3) 



f(p) even keeps an equilibrium shape after decoupling in 
two regimes || , namely ultrarelativistic decoupling (E ss 
pc, T — TjjRd/R, fi = hdRd I R), with 



f R (p) = {exp[(pc-n)/kT}±iy 



(4) 



and nonrelativistic decoupling (E — /i w p 2 /2m — jj,un> 
Mkin = M - mc2 = ^in.p{ R D / R) 2 , T = T D (R D /R) 2 ), so 



/at(p) = {exp [(//(2m) - Mkin)/fcT] ± l}~ 



(5) 



The distribution of fine-grained phase space density 
is conserved in time (Liouvilles theorem). This fact is 
expressed by Eq. (j|), and it means that the phase space 
distribution at decoupling can be directly related to mea- 
surements of dark matter phase densities today. Applied 
to conservation of the maximum phase space density this 
was the basis for the Tremaine-Gunn limit on dark mat- 
ter fermion masses 0] , later generalized to bosons (where 
no maximum exists) by means of the average phase space 
density §. 

Recently Hogan and Dalcanton Q reconsidered the is- 
sue of the primordial average phase space density com- 
pared with observations of dark matter phase space den- 
sities in halo cores. From the distribution function 
they calculated mass density p — mn, and pressure 
P = (g/h 3 )Jp 2 /(3E)fd 3 p « (g/h 3 )Jp 4 /(3mc 2 )4Trdp = 
mn < v 2 > /(3c 2 ), performing the calculation in the non- 
relativistic regime, where E «s mc 2 and p = mv. This 
leads to the following expression for the phase space den- 
sity, Q, defined by 



P 



U f{p)p 2 dpf' 2 

<v 2 >3 /2 h 3 71 [j /(p) p 4 dp ]3/2- [ > 

n S uses this expression in the 
= 1) 



Hogan and Dalcanto 
form (units with h — c 



lx = qxgxm x 



(7) 



for particle type X, where qx = 0.0019625 for a relativis- 
tically decoupling, u = fermion, and qx = 0.036335 for 
a relativistically decoupling, T = 0, /i » m degenerate 
fermion, where the g-values come from taking the com- 
plete integrals over the distribution function. 

For the remainder of this investigation we return to 
dimensional units and keep factors of h and c. In these 
units 



Qx 



] x h 3 m~ i g- 1 = 4tt 



[J f{p)p 2 dp} 5 / 2 
[J f(p) P 4 dp} 3 / 2 ' 



(8) 



For a degenerate fermion with relativistic decoupling 
(limit (m - n)/kT D — > -oo) this gives qRFdeg = 
4tt5 3 / 2 /3 5 / 2 = 9.0128415, so / = 9 /9.0128415 is in gen- 
eral a measure of the average occupation number in a 
phase space distribution (in units of g/h 3 ). For a zero 
chemical potential fermion the corresponding average q- 
values for relativistic and non-relativistic decoupling (in 
the latter case zero chemical potential means /ikm = 0) 
are quFO — 0.4868039 and qNFo = 1-9223, whereas the 
similar numbers for bosons are quso = 0.9071055 and 
qNBo = 21.521. Notice the larger values for bosons, that 
express the fact that bosons have a higher fraction of 
low momentum, high phase space density particles than 
fcrmions. 

So far only average characteristics of the distributions 
have been discussed, but it turns out to be quite inter- 
esting to study the whole distribution of phase densities. 
Figure 1 shows these distributions in a plot of q(p)/qx 
as a function of the fraction of particles with momentum 
less than p (calculated at decoupling, but the distribu- 
tion is conserved in time by Liouvilles theorem). Here 
qx is given in Eq. (Q) integrating from to oo, whereas 
q(p) is defined by the same equation, but integrating only 
from to p 1 10 1 . The ratio therefore illustrates the am- 
plification of phase space density relative to the average 
for a given dark matter distribution function if only the 
densest parts of phase space are utilized, for instance in 
formation of a galaxy halo. 

Notice that relativistically decoupling fermions show 
an order of magnitude amplification, whereas two orders 
of magnitude can be gained for relativistically decoupling 
bosons and several orders of magnitude for nonrelativis- 
tically decoupling bosons. 

Another way of illustrating the amplification effect is 
shown in Figure 2, where amplification is plotted as a 
function of the dimensionless momentum x = poc/kT^, 
which is the natural integration variable in the calcula- 
tions. For particles that are nonrelativistic at the epoch 
of galaxy formation, R g , x is related to particle speed v 
at that time by 



v mc 2 R g 
c kTu Rd 



(9) 



The occurrence of a phase space amplification fac- 
tor at low momenta is a natural consequence of 
Eq. (^). The fine-grained occupation number (and 
therefore also the coarse-grained phase space occupa- 
tion) has a maximum at p = equal to / ma x = 

{exp [(mc 2 — /Lt)/fcT] ± l} , which is 1 for degenerate 
fermions ((mc 2 — p,D)/kTu — > — oo), 0.5 for fermions de- 
coupling when (mc 2 — fiD)/kTD = 0, diverges for bosons 
in the same limit, and equals exp [(/io — mc^/feTo] for 
fermions and bosons in the limit (mc 2 — po)/kT]j — > oo. 
Quantitatively, the amplification factor for a relativis- 
tically decoupling fermion behaves like 9.26(1 — 5x/l6) 
to first order in x (Fig. §), or 9.26(1 - 0.691i^/ 3 ) ex- 
pressed in terms of the fraction of fermions, Fp (Fig. |l|). 
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The similar limits for bosons are 19.59x : (1 — 7x/30), 
■ 4.57 (notice that these factors diverge for 



or 8.93F, 



-1/2 



small x or Fb). 

Using entropy conservation in the cosmic expansion 
R g and Rd entering the equation for x are related by 
g* g T^ g Rg = g^oT^ D R 3 D , where g* counts the total num- 
ber of effective particle degrees of freedom at the given 
epoch. Today and at galaxy formation, g* g = 43/11. In- 
troducing the redshift of galaxy formation z g via T ig = 
(1 + z g )T 7 o, where the present photon temperature is 
T 7 q = 2.726K, x can be expressed as 
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(10) 



The value of g„D depends on the epoch of decou- 
pling. For standard neutrino decoupling at kTry ~ lMeV, 
g*D = 43/4. But much higher values of g^o are possible 
for earlier decoupling, with g^o w 50 above the quark- 
hadron phase transition temperature, kT « lOOMeV, 
g*D ~ 100 above kT rj 200GeV, and even higher val- 
ues possible earlier on. 

The value of not only determines the x-v rela- 
tion, but also crucially impacts on the total mass density 
contribution of the particle in question, and thereby de- 
termines its potential as a dark matter candidate. For 
the case of relativistic decoupling ((to — hd)/Td = 0, 
to/Tq — > 0) fermions contribute to the cosmic density 
like 



„2 \ 



n x h 2 = 0.0572 fel I g/g^D, 



(11) 



where h is now the Hubble-parameter in units of 100 
km s _1 Mpc -1 , and a similar expression (multiplied by 
4/3) applies for bosons. 

It is interesting to note, that an explanation of the 
highest phase space densities measured, those in dwarf 
spheroidals of order lO _4 M0pc -3 (km/s) -3 , for a g = 1 
boson requires a particle mass of only 224 eV (309 eV) if 
the densest 1% (10%) of the bosons are used, increasing 
to 681 eV for the average occupation number (masses are 
reduced by a factor 2 -1 / 4 if g = 2). Such masses are not 
in conflict with estimates of cosmic density for reasonably 
high values of g*D, and x can indeed be low enough for 
reasonable values of z g to select only the densest part of 
phase space for dwarf spheroidals with typical velocity 
dispersions below 10 km/s. 

For fermions with g — 2 the high-phase space density 
selection expected for formation of dwarf spheroidals re- 
duces warm dark matter particle mass limits from 669 eV 
to 383 eV, again loosening constraints on g^jj @- 

Even stronger effects may occur in a non-relativistic 
decoupling regime (upper dotted curve in Fig. 1), or one 
might consider adding extra features to the dark matter 
distribution functions, such as a small amount of Bose- 
Einstein condensation in the zero momentum, infinite 
phase density part of a boson distribution. 



Such "fine-tuning" may further loosen constraints on 
dissipationless dark matter, but even without it, part of 
the observed trend for the core phase space density to de- 
crease with increasing gravitational potential when going 
from dwarf spheroidal to dwarf spiral galaxies is naturally 
explained by the selection of low-momentum dark matter 
particles described here || ■ At the same time constraints 
on WDM particle masses and decoupling epochs are sig- 
nificantly reduced. These considerations should be taken 
into account and tested in detailed numerical simulations 
of dark matter halo formation, which are needed to set- 
tle the question of whether dissipationless particles may 
after all account for the dark matter in galaxies. 
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FIG. 2. As Figure 1, but as function of the dimensionless 
momentum, x. 



FIG. 1. Phase space density of dark matter particles 
in units of the mean density as a function of fraction of 
particles, F. Upper (lower) solid curves are fermions with 
(m — hd)/Td = and m/To — * (co). Lower (upper) dotted 
curves are for bosons in the same limits. Dashed curves are 
for fermions and bosons alike in the limits (m — hd ) /Td — * oo 
for m/To — > (upper) and m /Td — > oo (lower). Fully degen- 
erate fermions ((m — hd)/Td — > — oo) have no amplification 
factor, i.e. q/qx = 1. 
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